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Abstract 


A  theory  of  pump-probe  spectroscopy  is  developed  in  which  optical  fields 
drive  two-photon  Raman  transitions  between  ground  states  of  an  ensemble  of 
three-level  atoms.  Effects  related  to  the  recoil  the  atoms  undergo  as  a  result 
of  their  interactions  with  the  fields  are  fully  accounted  for  in  this  theory. 

The  linear  absorption  coefficient  of  a  weak  probe  field  in  the  presence  of  two 
pump  fields  of  arbitrary  strength  is  calculated.  For  subrecoil  cooled  atoms, 
the  spectrum  consists  of  eight  absorption  lines  and  eight  emission  lines.  In  the 
limit  that  xi  “C  X:i ,  where  xi  and  X'2  are  the  Rabi  frequencies  of  the  two  pump 
fields,  one  recovers  the  absorption  spectrum  for  a  probe  field  interacting  with 
an  effective  two-level  atom  in  the  presence  of  a  single  pump  field.  However 
when  xi  >  X2  >  new  interference  effects  arise  that  allow  one  to  selectively  turn 
on  and  off  some  of  these  recoil  induced  resonances. 

I.  INTRODUCTION 

Recent  advances  in  laser  cooling,  atom  optics,  and  Bose-Einstein  condensation  have 
underlined  the  role  played  by  atomic  recoil  in  atom-held  interactions.  A  measure  of  the 
importance  of  recoil  effects  is  the  recoil  frequency,  =  hk2/2M,  associated  with  the 
absorption,  emission  or  scattering  of  radiation  of  wavelength  A  =  27r//c  by  an  atom  of 
mass  M.  Once  this  quantity  becomes  greater  than  or  comparable  to  decay  rates  or  Doppler 
widths  that  characterize  the  spectral  response  of  atoms,  recoil  can  lead  to  new  features  in 
absorption  or  emission  line  shapes.  One  class  of  such  phenomena  has  been  termed  recoil- 
induced  resonances  (RIR)  [1-8],  which  occur  when  a  weak  probe  and  strong  pump  held 
simultaneously  drive  a  given  atomic  transition.  Interesting  in  their  own  right,  the  RIR  have 
been  used  to  determine  the  velocity  distribution  of  laser-cooled  atoms  f7j],  as  a  probe  of 
Bose-Einstein  condensates  [0]  flO]],  and  in  a  feedback  mechanism  in  stochastic  cooling  [111]]. 
Related  to  the  RIR  is  the  so-called  collective  atomic  recoil  laser  (CARL),  which  operates  on 
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similar  principles  but  in  a  somewhat  different  parameter  range  [|12j] .  Both  the  RIR  and  CARL 
represent  new  diagnostic  probes  of  cold-atom  sysytems.  Recently,  Moore  and  Meystre  f_I5J 
proposed  that  CARL  be  used  to  entangle  optical  and  matter  fields,  as  well  as  to  entangle 
different  modes  of  the  condensate  excited  by  the  optical  fields.  Our  discussion  is  limited 
to  situations  in  which  the  collective  effects  associated  with  CARL  can  be  neglected.  In 
this  paper,  we  combine  RIR  with  ground  state  spectroscopy  fid]  to  obtain  qualitatively  new 
features  in  the  probe  absorption  spectrum. 

The  scheme  we  adopt  is  based  on  the  model  developed  in  [|14j]  involving  Raman  transi¬ 
tions,  but  for  which  all  effects  associated  with  atomic  recoil  were  ignored.  In  that  work  a 
new  type  of  interference  was  discovered,  allowing  one  to  selectively  turn  on  and  off  certain 
lines  in  the  absorption-emission  spectrum  by  controlling  the  ratio  of  the  Rabi  frequencies  of 
the  two  fields  that  comprise  the  two-photon  pump  field.  Interference  in  a  dressed  state  basis 
occurs  between  pathways  involving  the  probe  field  and  each  of  the  two  pump  fields.  Since 
the  two  pump  fields  impart  different  recoil  momenta  to  the  atoms,  it  is  not  at  all  obvious 
that  the  interference  persists  when  recoil  splittings  are  resolved.  Part  of  the  motivation 
for  our  calculation  is  to  examine  this  question.  In  addition,  we  show  that  the  interference 
persists  even  if  the  pump  fields  are  in  quantized,  Fock  states. 

The  probe  absorption  spectrum  consists  of  as  many  as  eight  absorption-emission  doublets 
which  are  fully  resolvable  if  04,  >  7,  where  7  is  some  effective  ground  state  lifetime.  This 
is  in  contrast  to  the  RIR  spectrum  on  dipole  allowed  optical  transitions  [|I|],  where  at  most 
one  absorption  emission  doublet  is  resolvable  if  04,  <  7e,  where  7e  is  an  excited  state  decay 
rate.  One  might  question  the  need  to  increase  the  number  of  recoil  doublets  in  the  probe 
spectrum,  since  a  single  doublet  can  be  used  to  probe  recoil  effects.  We  show  that  the 
additional  recoil  structure  reflects  interesting  quantum  dynamics  of  the  combined  atom-field 
system,  as  well  as  providing  some  new  applications. 

In  Sect.  II,  a  model  is  developed  for  the  interaction  of  the  atoms  with  the  pump  fields; 
dressed  states  of  the  atom  plus  pump  fields  are  defined.  In  Sect.  Ill,  the  interaction  with  the 
probe  field  is  introduced  and  the  dressed  state  picture  is  used  to  obtain  the  probe  absorption 
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spectrum  in  the  secular  limit.  In  Sect.  IV  we  discuss  the  results  and  possible  applications. 
Nonsecular  contributions  to  the  absorption  coefficient  are  calculated  in  an  Appendix. 

II.  EQUATIONS  OF  MOTION 

The  atom-field  system  is  indicated  schematically  in  Fig.l.  Ground  state  levels  |1  >  and 
1 2  >  are  pumped  incoherently  with  rates  Ai(p)  and  A2(p),  respectively,  and  the  population 
of  both  states  decay  with  rate  7.  If  states  |1  >  and  |2  >  represent  stable  ground  states  of 
the  atom,  then  the  pumping  rates  and  decay  rate  constitute  a  simple  model  for  atoms  that 
enter  and  leave  the  interaction  volume.  The  ground  to  excited  state  transition  frequencies 
are  denoted  by  uej  (j  =  1,2).  The  pump  fields  1  and  2,  which  constitute  the  two-photon 
pump  field,  are  denoted  by  the  coupling  constants  57  and  r/2,  respectively.  Pump  field  1 
couples  only  state  |1  >  and  excited  state  |e  >,  while  pump  2  couples  only  states  \2  > 
and  \e  >.  The  pump  fields  have  frequencies  Qi  and  02,  and  propagation  vectors  ki  and 
k2,  respectively.  In  this  section,  equations  are  derived  for  the  atom-pump  field  interaction, 
neglecting  the  incoherent  pumping  and  decay  of  the  ground  state  levels.  In  the  following 
section,  the  interaction  of  the  atoms  with  the  probe  field,  Ep,  which  couples  only  states  |1  > 
and  \e  >,  is  incorporated  into  the  calculation,  as  are  the  incoherent  pumping  rates  Ai(p) 
and  A2(p),  and  the  ground  state  decay  rate  7. 

In  contrast  to  [114]  but  as  in  [l3j] ,  we  use  a  quantized  description  of  the  pump  fields.  If  the 
pump  field  detunings, 


A2  —  Q2  —  ue2  ~  Ai  —  Qi  —  u>ei  =  A,  (2) 

are  sufficiently  large  such  that  |xt,2/ A|2  -C  1  and  1 7e / A |  -C  1,  where  the  77,2  are  defined 
below  and  qe  is  the  excited  state  decay  rate,  it  is  possible  to  adiabatically  eliminate  the 
excited  state  to  arrive  at  an  effective  Hamiltonian  involving  only  states  |1  >  and  \2  >  which 
is  of  the  form 
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H  —  Ha  +  Hr  +  Har ; 

„2 


p 


(fiiJi  H - )  1 1,  p  ><  1,  p|  +  (hiL>2  H - )|2,  p  ><  2,  p| 

2  M  2  M 


Hr  =  hfiia\ai  +  M22c4a2; 


t. 


Har  =  YJ 


h 


M*eikl2'R 


2,  p  ><  1,  p|c4ai  +  ^ 


g2g*e^n 


1,P  ><  2,p|a|a2 


p 


(3a) 

(3b) 

(3c) 

(3d) 


A  1  1  "  A 

where  Ha  is  the  Hamiltonian  for  the  atom  in  which  the  center-of-mass  momentum  p  has 
been  quantized  using  periodic  boundary  conditions  in  a  volume  V  (assuming  that  the  atoms 
are  free  and  not  subject  to  some  trapping  potential),  Hr  is  the  free  held  Hamiltonian  for  the 
two  pump  fields,  Har  represents  the  interaction  of  the  two-photon  pump  held  with  an  atom, 
in  the  rotating-wave  approximation,  and 

kp  =  k,  k  j. 

Bare  states,  \j,  p;  ni,  n2  >,  are  defined  to  be  eigenstates  of  Ha  +  Hr,  where  j  =  1,  2  labels  the 
internal  state  of  the  atom  and  n i  and  rt2  are  the  number  of  photons  in  pump  helds  1  and  2, 
respectively.  A  term  in  the  Hamiltonian  corresponding  to  the  light  shifts  of  the  ground  state 


levels,  His  =  P  ><  1,  pK^i  +  p  ><  2,  p|4a2 

Eq.  (Saj);  such  light  shifts  can  be  included  by  a  redefinition  of  the  ground  state  frequencies, 


\92\2 


,  has  been  omitted  in 


|92|2 


,t. 


Isil2 


,t. 


cu2  d-  2n2  J  — *■  u>2  ana  u i  +  <^a[a i )  — >  cui. 

The  matrix  elements  of  the  operator  e*k  R  in  the  momentum-state  basis  are 

<  p|ejk'R|p'  >=<  p|p'  +  Kk  >=  hp,p'+nk- 


This  allows  one  to  rewrite  the  interaction  term  as 


Har  =  YJ 


n^-|2,p  +  hki2  ><  l,p|4ai  +  1 1, p  ><  2,  p  +  hki2|a|a2 


.929! 


p 


(4) 


(6) 


The  Hamiltonian  H  results  in  an  inhnite  ladder  of  decoupled  two  state  manifolds  (p,  ni,  n2) 
involving  the  states  |l,p;ni,n2  >  and  |2,p  +  hki2;ni  —  l,n2  +  1  >.  The  Hamiltonian  for 
the  manifold  (p, ni,n2)  is 


H{  p,  ni,n2)  =  e(p,  m,  n2)I  +  h 


1  -S(  p)/2  G*  ^ 
K  G  5(p)/2 } 


(7) 
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where  I  is  the  identity  matrix 


G  = 


=  \G\ei<t>d, 


e(p,ni,n2)  =  ft(nifii  +  n202)  +  ^(<^p  +  <np+mci2  +  cvi  +  cn2  +  fl2  -  Qi),  (8) 


pk 
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<5(p)  -  hr2  -  uhkl2  M  , 

<5i2  =  Ai  —  A2  =  (Jli  —  S72)  —  o;2i, 


=  S’  and  *2  =  ^2\/n2  +  1,  xi  =  9\ \fn\- 

The  dressed  states  are  dehned  to  be  the  eigenstates  of  the  matrix  in  Eq.(!7|),  with  energies 

j?  ,  \  i  \  I  ^ab( p) 

Eb,a( P)  =  £(P ,ni,n2)  ± - - - ; 


uab(  p)  =  v/  4|G|2  +  <5(p)2, 


and  associated  eigenkets  [Il5|] 


where 


and 


V 


I -Bn  > 


=  T*(p) 


^  |2,  p  +  hki2;  m  -  1,  n2  +  1  >  ^ 


V 


|l,p;ni,n2  > 
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T(P)  = 


7  ei<t>d/ 2  cos  —  e  Ad/2  sjn  ^ 


eAd/2  sin  [5)(p)]  e  Ad/2  cog 


cos  [0(p)]  = 


1  1  +  Ap) 

2  l  oub(p) 


1/2 


(9) 


(10) 


(12) 


The  value  of  0(p)  is  restricted  to  the  range  0  <  0(p)  <  7t/4  for  <5(p)  >  0  and  7t/4  <  0(p)  < 
7t/2  for  <5(p)  <  0.  For  0(p)  ~  0  (<5(p)  >  0  and  <5(p)/|G|  S>  1),  |A0  >~  1 2,  p  +  hki2;  n\  —  1,  n2  + 
1  >  while  for  0(p)  ~  7t/2  (5(p)  <  0  and  <5(p)  /\G\  1),  |B0  >~  1 2,  p+hki2;  m  —  1,  n2  +  l  >. 
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III.  PROBE  FIELD  ABSORPTION  IN  THE  SECULAR  LIMIT 


The  effect  of  the  probe  held  is  to  induce  transitions  between  states  in  different  manifolds. 
As  is  customary  in  dressed  atom  approaches,  the  probe  is  treated  as  a  classical  held, 

E(R,t)  =  -?Epel(k‘’'R-a^  +  c.c. ,  (13) 

A 

where  'e  is  a  unit  polarization  vector.  For  our  problem,  however,  this  choice  represents 
a  hybrid  approach,  since  two-quantum  processes  involving  the  probe  held  and  either  of  the 
pump  helds  mix  classical  and  quantized  helds.  Although  the  probe  held  is  treated  classically, 
its  effect  on  the  momentum  of  the  states  must  be  accounted  for  explicitly.  If  the  probe  held 
detuning  on  the  1  — >  e  transition  is  sufficiently  large  to  be  consistent  with  the  adiabatic 
elimination  of  the  excited  state,  all  transitions  involving  the  probe  occur  via  two-quantum 
transitions  involving  the  probe  held  and  either  of  the  pump  helds.  In  the  bare  state  basis, 
starting  from  the  (p,  711,712)  manifold,  probe  held  absorption  corresponds  to  transitions 
|l,p;7i!,7i2  >— >  1 1,  p  +  hkp^ri!  +  1,77.2  >  or  |l,p;ni,n2  >— >  |2,  p  +  Tikp2\  nun2  +  1  >, 
where  the  second  photon  is  emitted  into  either  the  pump  1  or  pump  2  modes,  respectively. 
Similarly,  probe  gain  corresponds  to  transitions  1 1,  p;  711,  n2  >— >  1 1,  p  —  hkpi;ni  —  1,  n2  >  or 
1 2,  p  +  hki2;  7ii  —  1,  ti2  +  1  >— >  1 1,  p  —  hkpi,  711  —  1,  n2  >.  This  is  illustrated  in  Fig.  2(a). 

This  picture  of  probe  held  absorption  or  emission  allows  us  to  reach  an  important  con¬ 
clusion  concerning  interference  between  pathways  involving  both  pump  helds.  The  two 
absorption  processes  shown  in  Fig.  2(a)  involve  different  hnal  states  and  do  not  interfere, 
reinforcing  the  possibility  mentioned  in  the  Introduction  that  interference  may  be  suppressed 
when  recoil  is  taken  into  account.  However,  both  hnal  states  belong  to  the  same  manifold  - 
the  (p  +  hkpi,  711  + 1,  ti2)  manifold.  As  such,  when  these  states  are  dressed  by  the  two-photon 
pump  held,  each  state  in  the  hnal  state  manifold  will  be  coupled  to  the  initial  state  by  two 
separate  pathways  involving  the  probe  held  and  each  of  the  pump  helds.  This  implies  that 
terms  in  the  probe  absorption  depending  on  the  simultaneous  presence  of  both  pump  helds 
will  exhibit  interference  effects.  In  the  case  of  probe  gain,  it  is  the  initial  states  that  differ, 
but  the  overall  conclusion  remains  unchanged. 
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The  dressed  state  approach  provides  a  convenient  and  relatively  easy  method  for  obtain¬ 
ing  the  probe  absorption  spectrum  in  the  secular  limit,  where  the  frequency  separation  of 
the  dressed  states  in  a  given  doublet  is  much  larger  than  the  ground  state  decay  rate  7.  [It 


is  assumed  from  this  point  onward  that  ye  (xi)2/A)2  <C  7,  implying  that  the  dressed  states 


decay  with  rate  7  [l3j].  One  need  calculate  only  the  dressed  state  energies  and  the  transition 


matrix  elements  to  obtain  the  spectrum.  A  more  detailed  treatment  of  the  problem,  allowing 


one  to  calculate  non-secular  contributions,  is  presented  in  Appendix  B. 

It  is  straightforward  to  generalize  the  dressed  states  defined  in  the  previous  section 
to  include  the  two  manifolds  coupled  to  the  initial  manifold  by  the  probe.  The  0, 
1,  and  2  manifolds  refer  to  (p,  77,1,71,2)  =  { 1 1,  p;  77,1,  7t,2  >,  |2,  p  +  hki2]  rq  —  1, 77,2  +  1  >}, 
(p  +  hkpi,?7i  +  1,  ?t,2 )  =  {|l,p  +  Kkpi]rii  +  l,n2  >,  |2,p  +  ^,(ki2  +  kpi);77i,n2  +  1  >}  ,  and 
(p  -  hkpi,  77,1  -  l,ra2)  =  { 1 1,  P  -  hkpl;ni  -  1  ,n2  >,  |2,p  +  h( k12  -  kpl);ni  -  2,n2  +  1  >}, 
respectively.  Taking  the  central  energy  of  the  initial  manifold  (p,  77,1,  n2)  arbitrarily  equal  to 
zero,  one  finds  that  the  dressed  state  energies  are  given  by 


(14) 


where 


^ab(p)  =  vW+W, 


(15) 


(p  +  hkpi)  ■  k12 

M 

(p  -  hkpi)  ■  kj2 

M 


(16) 


h2(p)  =  h(p  -  hkpi)  =  S 1 2  -  07Sk12 
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and  it  has  also  been  assumed  that  n\,n2  1  such  that  G ^  =  g*g2 a/ui  (n2  +  1)/A  ~ 

G(1)  =  g*g2\/Jn\  +  1)  {n2  +  1)/A  m  G(2)  =  g*g2\/Tji\  -  1)  (n2  +  1)/A  =  G.  The  dressed 
state  angles  are  given  by 


cos  [0i(p)] 


1 

2 


(19) 


and  dressed  state  kets  are  dehned  by 


\Ai  > 
I  Bi  > 


\ 
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/ 


=  T*(p) 


2,  P  +  h  [ki2— (— l)1  (1  -  <5i,o)  kpi] ;  ni  -  1  -  (-1)*  (1  -  ^i0) ,  n2  +  1  > 
1 1,  p— (— 1)*  (1  -  <5i,0)  ftk.pi,  n\  -  (-1)*  (1  -  <5j,0) ,  n2  > 


\ 


(M) 


where 


g^d/2  cos  [0j(p)]  — e  sin  [0*(p)] 

^  e^d/2  sin  [0j(p)]  e~l<l>d/2  cos  [0j(p)]  j 


(22) 


The  absorption  coefficient,  a ,  is  proportional  to  the  rate  at  which  energy  is  absorbed 
from  ( a  >  0)  or  emitted  into  ( a  <  0)  the  probe  held.  Absorption  corresponds  to  transitions 
from  initial  dressed  states  /  =  Aq,  B0  to  final  state  dressed  states  J  =  A\,B\,  while  emis¬ 
sion  corresponds  to  transitions  from  initial  dressed  states  /  =  A0,  B0  to  final  state  dressed 
states  J  =  A2lB2.  For  a  given  transition,  the  contribution  to  the  absorption  coefficient  is 
proportional  to  NVhQp'ypjj /\Ep\\  where  N  is  the  atom  density  and  pjj  is  the  steady  state 
population  in  state  J  owing  to  the  /  — *  J  transition.  The  final  state  population  pjj(p)  is 
equal  to  [Aj(p)/7]  |  <  J\VP\I  >  \2Ljj,  where  A/(p)  is  the  pumping  rate  for  initial  dressed 
state  /  [|I6|],  <  >  is  a  matrix  element  for  the  I  —  J  transition  and  L/j  is  a  Lorentzian 

having  width  7,  centered  at  the  I-J  transition  frequency.  The  transition  frequencies  A/j(p) 
may  be  directly  read  from  Fig.  2(b)  or  obtained  from  Eq.  (|14j) .  Consequently,  it  is  neces¬ 
sary  to  specify  A/(p)  and  to  calculate  the  transition  matrix  elements  <  >  in  order 

to  arrive  at  an  expression  for  the  probe  absorption. 

In  the  bare  state  representation,  the  pumping  matrix  for  the  intitial  state  manifold  is 
taken  to  be  of  the  form 
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A(P)  = 


^Pn(p) 


0 


V 


o  p22(p+^k 


12 


\ 


) 


pump 


Ai(p)  0 
0  A2(p  +  ^ki2 


\ 


which,  when  converted  to  the  dressed  state  basis  becomes 

Ad(p)  = 


^  Aj4(p)  Aab(p)  ^ 


Aab(p)  Ab(p) 


J 


A2(p  +  hk12)  cos2  [0o(p)]  +  Ai(p)  sin2  [0o(p)]  \  [A2(p  +  Kkr2)  -  Ai(p)]  sin  [26>0 (p)] 

^  \  [A2(p  +  hki2)  -  Ai(p)]  sin  [200 (p)]  Ai(p)  cos2  [0o(p)]  +  A2(p  +  hk12)  sin2  [0o(p)] 

In  the  secular  limit,  the  off-diagonal  pumping  terms  can  be  neglected  since  they  give  rise  to 
terms  of  order  7/0;^ (p)  <C  1.  Moreover,  for  the  present  we  will  set  A2(p  +  fik\2)  =0  and 
generalize  the  results  to  nonvanishing  A2(p  +  hki2)  in  the  next  section.  Thus,  we  take  a 
pumping  matrix  of  the  form 


/ 


(p)  = 


\ 


V 


Ai  (p)sin2  [0O  (p)]  0 

0  Ai  (p)  cos2  [0O  (p)]  J 

We  now  turn  our  attention  to  the  transition  matrix  element  <  J\VP\I  >  .  As  an  ex¬ 
ample  consider  <  Ai|h^|A0  >,  which  is  the  amplitude  for  the  transition  \A0  >— >  \A\  > 
involving  probe  absorption  and  pump  1  or  pump  2  emission.  This  transition  is  illustrated 
in  Fig.  2(c).  The  probe  couples  only  to  the  |l,p;rai,n2  >  part  of  \A0  >,  leading  to  a  factor 
_e*9W2  sin  [0o(p)]_  The  absorption  of  the  probe  is  followed  by  emission  into  pump  2  taking 
the  atom  to  the  1 2,  p  +  h( ki2  +  kpl );  ri\ ,  n2  +  1  >  component  of  |  A\  >  and  leading  to  a  factor 
g-^d/2  cog  [01(p)j  or  emission  into  pump  1  taking  the  atom  to  the  the  1 1,  p  +  hkpi;  n\  +  l,  n2  > 
component  and  leading  to  a  factor  —el^d/2  sin  [0i  (p)] .  The  coupling  strengths  for  these  two, 
two-photon  transitions  are  G\  and  G 2,  respectively,  where 


G 1  = 


XpX  1  ^  XpX  2 

Cj2  — 


A  ’ 


(25) 


The  two  processes  add  coherently  and  one  finds 


<  A^VIAq  >=  h  {G*2e^d  cos  [0i(p)]  -  G{  sin  [0i(p)]}  {-  sin  [0o(p)]}  e~iQA. 
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Other  matrix  elements  are  calculated  in  a  similar  manner.  For  probe  gain,  pump  fields  1  and 
2  couple  to  the  1 1,  p;  ni,  n2  >  and  |2,  p+hki2;  n\  —  1,  n2  +  l  >  components  of  the  dressed  states 
in  the  0  manifold,  respectively,  while  the  probe  field  couples  to  the  1 1,  p  —  hkp i ;  ri\  —  1,  n2  > 
component  of  the  dressed  states  in  the  2  manifold.  Explicit  expressions  for  the  matrix 
elements  are  given  in  Appendix  A. 

Combining  all  transitions  and  summing  over  p,  one  finds  an  absorption  coefficient  pro¬ 
portional  to 


E  E 


Mp) 


I={A0,B0} 


J={A\,B\} 


<  J\VP\I  >  |2 

-  A/j(p)]2  +  72 


E 

J={A2,B2} 


\<j\vp\i>\ 2  \ 

[ap  -  A/j(p)]2  +  7 2  ) 

(26) 


For  a  sub-recoiled  cooled  vapor,  we  can  set  Ai(p)  =  Ai5Pi0,  such  that  A^(p)  =  A^tfp^o  and 
Ab(p)  =  AS(5Pi0,  where 


A  a  =  Ai  sin2  60]  As  =  Ai  cos2  90 , 

6i  =  9i{ 0)  [and,  for  future  reference,  A jj  =  A/j(0),  c o^B  =  u>^B(0),  etc.]  and  5Pio  is  a 
kronecker  delta.  In  this  limit,  one  finds  the  absorption  coefficient  in  the  secular  limit  to  be 
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LaqA-l  (A7) 
LsoAi  (A') 

AaoBi(A') 

^BoSi(A') 

La0a2  (A  ) 
LaqB2  (A7) 
Lb0a2( A') 
Lb0b2( A') 

where 


H 

I A  | 


+ 


^00 sin  #i  -  ^  cos  ^  sin2  OqLaoA!  (A')  +  ^  cos2  OqLbqA!  (A')^ 

V>0  cos  0i  +  -  sin  0^  f —  sin2  90LAob1  (A')  +  —  cos2  90LBob1  (A7)^ 

V  J  \  7  7  ) 


ijjr]  sin  0O - cos  0o 

V 


—  '0r/  cos  0o  H —  sin  0O 

V 

72 


A /i 

7 

As 

7 


(sin2  92LAoa2( A')  +  cos2  92LAob2(A')) 
(sin2  02LBoj42(A/)  +  cos2  92LBob2(A')) 


(A'  -  ^kpl  -  i(^!S  -  w^))  +  72 

_ 7^ _ . 

(A'  -  ^kpl  +  |  (^  +  w^b))  +  72 

_ 7^ _ . 

(A'  -  ^kpl  -  |  (o>2b  +  W AB ))  +  72 

_ 7^ _ . 

(A'  -  ujhkpl  +  \  (ufB  -  w^b))  +  72 
2 
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(A/  +  ^kpl  +  |  (w2b  -  w2b))  +  72 

_ 7^ _ . 

(A7  +  ^kpi  + §  (<^i°i  +  +  72 

_ 7^ _ . 

(A7  +  ^hkpl  ~  \ (^ab  +  wab))  +  72 
72 


aii,,  i  (. ,(°)  _  ,  ,(2) 

^  +  k-7kpi  2 


(27a) 

(27b) 

(27c) 

(27d) 

(27e) 

(27f) 

(27g) 

(27h) 

(27i) 


<l>  =  A/|A| 

is  the  sign  of  the  detuning  for  each  of  the  three  fields, 

V  =  \/ 1  Xi  I  / 1 X2 1 , 


12 


A  /  A  ^'kp  1  •  ki2 

A  -  ~  2M  ’ 

8p\  D p  Hi,  3p2  D p  D2  ^21 ) 


«o  — 


kpNd\e 

2he0j  ’ 


and  die  is  a  bare  state  dipole  moment  matrix  element  for  the  |1)  — >  |e)  transition. 


The  hrst  four  resonances  in  (  — 

V  «o 


correspond  to  probe  absorption  while  the  last  four 


correspond  to  probe  gain.  The  spectrum  is  shown  in  Fig.  3.  The  line  widths  of  all  the 
resonances  equal  7;  consequently,  the  recoil  induced  resonances  in  ground  state  spectroscopy 
are  fully  resolved  if  07  >  7.  The  secular  contribution  does  not  vanish  in  the  limit  that 
|<5(p  =  0)/Gj  <C  1,  even  though  =  A B  in  this  limit.  As  long  as  the  recoil  frequency  is 
larger  than  7  and  the  Doppler  width  associated  with  the  two-photon  pump  transition,  the 
absorption  and  emission  contributions  to  the  probe  response  do  not  cancel  one  another. 


The  most  significant  feature  of 


is  that  the  line  strengths  involve  factors  such  as 


^iprj  sin  0\  ~  2  cos  which  allows  one  to  manipulate  the  strength  of  the  lines  by  controlling 
the  sign  of  the  held  detuning  and  the  ratio  of  the  pump  held  amplitudes.  These  factors  are  an 
indication  of  interference  between  the  two  ”  two-photon  probe”  helds  which  can  both  lead  to 
absorption  or  gain  in  Ep.  Because  the  two  lines  in  the  doublets  have  different  strengths,  one 
can  adjust  ip  and  77  to  turn  oh  one  of  the  lines.  For  example,  the  absorption  doublet,  La0a1 
and  LBqBi,  consists  of  the  lines  at  A'  =  ujhkpl  +  K^ab~ujab)  and  ^Kkpl  -  K^ab  ~ 00 ab)  with 
strengths  ~  ^77  sin  6\  —  ±  cos  0i^  and  ~  ^77  cos  Q\  +  ^  sin  9\  j  ,  respectively.  Choosing 
7 p  =  +1  and  r/2  =  cot  turns  ”oh”  the  hrst  line  while  ip  =  —1  and  r]2  =  tan  9\  turns  ”oh” 
the  second  line.  This  is  shown  in  Fig.  4.  When  |Gj  and  812  are  much  larger  than  any  of 
the  recoil  terms,  0O  ~  0i  ~  02  and  the  emission  lines  are  also  turned  ”oh”.  Consequently, 
by  choosing  ip  =  +1  and  rj2  =  cot  6\  to  turn  oh  the  La0a1  and  LBqj 41  absorption  lines,  the 
La0a2  and  La0b2  emission  lines  are  also  turned  oh. 

A  particularly  interesting  case  occurs  when  ki  ~  k2  so  that  c =  u^B  =  ^ab  and  0o  = 
0i  =  02-  This  would  correspond  to  a  two-photon  pump  held  which  imparts  no  momentum 
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to  the  atoms  so  that  the  recoil  splitting  in  the  absorption  spectrum  can  be  attributed  solely 
to  the  recoil  due  to  the  probe  field  acting  with  either  of  the  pump  fields,  kp2  =  kp] .  In 
this  case,  the  line  La0a1  (A')  is  degenerate  with  LbqBx  (A')  and  La0a2  (A')  is  degenerate 
with  Lb0b2  (A')  .  Consequently,  the  spectrum  consists  of  three  absorption-emission  doublets 
centered  at  A'  =  0,  +u>abi  ~ujab-  Moreover,  the  lines  within  each  doublet  are  split  by  2ujfikpl 
which  is  independent  of  the  strength  or  detuning  of  the  pump  fields. 

When  the  effects  of  atomic  recoil  are  neglected  by  setting  all  recoil  momenta  to  zero  in 
one  obtains  the  same  absorption  spectrum  given  in  pf4j].  In  the  limit  that  G\  =  0,  one 
recovers  a  simple,  recoil  shifted  Raman  spectrum.  In  the  limit  that  G2  =  0,  one  recovers 
the  central,  secular  components  of  the  pump-probe  spectrum  associated  with  a  single ,  two- 
level  optical  transition  [J3j].  In  the  limit  that  77  <C  1,  while  G/ 7  1  remains  constant,  the 

absorption  spectrum  mirrors  that  for  the  pump-probe  spectrum  associated  with  a  two-level 
optical  transition  [13]. 


IV.  DISCUSSION 

For  a  subrecoil  cooled  vapor,  pumping  to  state  2  at  a  rate  A2  (p  +  hki2)  =  A2(5p+^ki2,o 
doubles  the  number  of  absorption  and  emission  lines  in  the  probe  spectrum,  but  does  not 
result  in  any  qualitatively  new  features.  There  will  be  an  additional  contribution  to  Eq. 
(J>2!)  in  which  9i  =  0*(O)  is  replaced  by  9t  =  9i(-hk12),  Au  by  Au  =  A7J(-hki2),  ulAB 
by  u j1ab  =  uab(~ hki2),AJ4  by  A2cos2$o  and  A b  by  A2sin20o-  The  absorption  coefficient 
contains  sixteen  lines  in  all.  The  eight  new  spectral  components  display  the  same  properties 
as  the  original  eight  but  are  displaced  by  an  amount  Figure  5  shows  the  secular 

absorption  spectrum  with  all  sixteen  components  when  A2  =  Aj.  In  the  absence  of  recoil, 
the  secular  absorption  coefficient  vanishes  when  A2  =  Ai  1114].  However,  when  recoil  is 
included,  the  absorption  and  emission  contributions  to  the  probe  response  do  not  cancel  one 
another  when  A2  =  Ai,  provided  the  recoil  frequency  is  larger  than  7  and  the  Doppler  width 
associated  with  the  two-photon  pump  transition. 
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In  Ref.  [jl4j]  interference  effects  similar  to  those  in  Eq.  (27a|)  were  found.  We  have  shown 
that  the  interference  is  independent  of  the  field  statistics  and  persisits  even  when  recoil 
induced  resonances  are  resolved. 

The  RIR  offer  several  possibilities  for  applications.  The  existence  of  a  number  of  tunable, 
well-resolved,  gain  peaks  allows  one  to  envision  experiments  in  which  lasing  occurs  at  one 
or  more  of  the  gain  positions.  These  peaks  could  be  adjusted  to  coincide  with  modes  of  a 
ring  cavity,  for  example.  By  sweeping  the  ratio  of  the  two  pump  field  intensities,  one  has 
a  mechanism  for  modifying  the  probe  gain  or  absorption,  to  the  point  of  total  suppression. 
The  central  frequency  of  the  absorption-emission  doublets  can  be  controlled  via  pump  field 
strength  and  detuning.  For  parallel  pump  fields,  the  frequency  separation  of  the  absorption 
emission  doublets  can  be  as  large  as  804,  and  is  independent  of  pump  field  strength  and 
detuning.  With  line  widths  approaching  1  Hz  or  less,  the  probe  spectrum  can  be  used  to 
measure  the  recoil  frequency  to  a  precision  of  order  ICE7;  this  precision  can  be  increased  if 
the  two-photon  pump  held  is  replaced  by  a  pair  of  counterpropagating,  two-photon  pump 
fields  [117] .  The  narrow  resonances  can  also  be  used  in  schemes  for  obtaining  ’’slow  light” 

PS- 

Pump-probe  spectroscopy  of  Bose-Einstein  condensates  represents  an  interesting  appli¬ 
cation  of  the  ideas  presented  here.  Bragg  spectroscopy  has  recently  been  demonstrated  in 
condensates  [H]  [|10l]  as  well  as  the  stimulated  generation  of  matter  waves  in  a  condensate 
by  Rayleigh  scattering  [119] .  Pump-probe  spectroscopy  using  electronic  excited  states  would 
be  unfeasible  in  condensates  since  there  are  no  stable  trapped  condensates  with  electronic 
states  which  may  be  populated.  Currently,  the  only  multi- component  condensates  consist  of 
two  hyperhne  states  in  8‘  Rb  [121]  and  the  Zeeman  states  in  the  F  —  1  manifold  of  optically 
trapped  Na  [122] .  Consequently,  pump-probe  spectroscopy  would  necessarily  involve  Raman 
transitions  between  stable  ground  states  in  the  manner  proposed  here.  The  RIR  spectrum 
of  a  weakly  interacting  Bose  condensate  should  yield  information  about  the  spectrum  of 
elementary  excitations  in  the  condensate  which  for  small  momenta  have  a  linear  dispersion 
relation  while,  for  large  momenta,  have  a  quadratic  dispersion  similar  to  that  of  free  atoms, 
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but  with  a  shift  due  to  the  mean-field  interactions  in  the  condensate.  In  addition,  the  line- 
widths  of  the  RIR  spectrum  should  be  given  by  the  zero  point  motion  of  the  condensate  in 
the  trapping  potential  provided  (yyyy)  kP\,2  >  7  and  Ax  is  the  size  of  the  condensate.  How¬ 
ever,  the  direct  application  of  the  results  presented  here  to  a  condensate  would  be  erroneous 
since  a  correct  calculation  of  the  RIR  spectrum  would  have  to  account  for  the  mean-field 
interactions  between  the  atoms.  This  will  be  pursued  in  future  work. 
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VI.  APPENDIX  A  -  MATRIX  ELEMENTS  OF  VID 

The  interaction  with  the  probe  field  may  be  expressed  in  terms  of  an  effective  two-photon 
interaction  Hamiltonian  similar  to  that  in  Eq.  (!3cl;), 

V/  =  f^^-e~^pt  |2,p  +  hkp2)  (l,p|  a\  +  |l,p  +  hkpi)  <1,  p|  a[  +  h.c^j 

(28) 

where  \p  —  yy deiEp  is  the  Rabi  frequency  for  the  probe  field  and  de i  =  <  e|d-e|l  >  is  a 
dipole  matrix  element.  The  matrix  representation  of  V/  with  respect  to  the  bare  state  basis 
has  the  following  nonvanishing  elements 

(l,p;ni,n2|  V/  |2,p  +  hk12  +  hkpl]  ni,  n2  +  1)  =  hG2elQpt] 

(l,p;ni,n2|  V/  |l,p  +  hkp^ni  +  1  ,n2)  =  hGieinpt ;  (29) 

(l,p-/fkpi;ni  -  l,n2|  V/  |2,p  +  hki2;ni  -  l,n2  +  1)  =  hG2elQpt; 

(l,p— ftkpi;m  -  l,n2|  V/  1 1,  p;  ni,n2)  =  hGieinpt ; 
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and  the  hermitian  conjugates  of  Eqs.  (29).  In  Eqs.  (29),  G\  and  G2  are  two-photon  probe 
Rabi  frequencies  dehned  as  G\  =  x*xi/A;  G2  =  x*PX 2/A.  Note  that  the  coupling  of  the  1 
and  2  manifolds  to  manifolds  other  than  the  0  manifold  has  been  ignored. 

The  matrix  V/d  represents  the  interaction  with  the  probe  held  in  the  dressed  state  basis 


and  is  dehned  as 

VId  =  TV;Tt;  (30) 

where  T  is  given  by  the  block  diagonal  matrix 

T  =  diag  [Ti(p),  T0(p),  T2(p)] .  (31) 

The  Rabi  frequencies  xP->  Xu  and  X2  may  be  expressed  as 

XP  =  \xP\e^  Xi  =  IXile^1;  X2  =  |X2|e*02;  (32) 

so  that  4>d  =  4>2  ~  4>i  +  f  (1  —  0)  since  G  =  \G\el^d.  The  matrix  elements  are: 

<  AilV/14)  >=  ^(^l)e-^sin0o(-|G'2|cos0i  +0|G'i|sin01);  (33a) 

<  Ri|V/|Ro  >=  e~iVlpt  cos  ( | <^2 1  cos  9\  —  i/j\Gi\  sin  ffi);  (33b) 

<  Bi\V  i\A0  >=  hei^~(pl)e-iQptsm90(-'ilj\G1\cosd1  -  |G2|sin0i);  (33c) 

<  Ri|V/|R0  >=  e~lQ‘pt  cos0o(0|Ca|  cosffi  +  |G2|sinffi);  (33d) 

<  Ro|V/|R2  >=  sin 02(-\G2\  cos  90  +  0|Gi|  sin0o);  (33e) 

<  R0|V/|R2  >=  e~iQpt  cos  02(|Cr2|  cos  90  —  0|Gh|  sin  d0);  (33f) 

<  R0|V/|R2  >=  he^-^e-iQpt  sin  92(-il;\G1\ cos  90  -  |G2|sin0o);  (33g) 

<  R0|V/|R2  >=  e-'lQ,pt  cos02(0|Gi|  cos  90  +  |G2|sind0)-  (33h) 

The  other  elements  follow  from  the  hermiticity  of  V 


VII.  APPENDIX  B 

In  this  Appendix  the  absorption  coefficient  for  the  probe  held  is  calculated,  without 
making  the  secular  approximation.  The  absorption  coefficient,  a,  and  index  change,  An, 
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arise  from  the  imaginary  and  real  parts  of  the  macroscopic  polarization  in  the  Maxwell-Bloch 
equations  for  the  probe  held.  They  are  given  by  the  expressions 


a  =  «o  Im 


7 Pie 

V* 
A  n 


V 


An  =  — «o  kp  1  Re 


7  Pie 

V* 

An 


V 


«o  — 


kpNdie 

2^eo7 


(34a) 

(34b) 

(34c) 


where  V  is  the  volume  and  p'le  is  the  part  of  the  bare  state  density  matrix  element  pie(R,  t) 
which  is  proportional  to  e-7k,,-R-n?Jp^  which  we  denote  by  //le(R,  £). 

Before  proceeding,  we  note  that  in  this  Appendix  all  summations  over  momentum  states 
have  been  converted  to  an  integration  over  a  continuum  of  states  via  the  standard  substi¬ 
tution  — >  A  f  d3p.  The  coefficient  p\e  is  related  to  the  coherence  in  position  space, 

p'le(R,  t),  and  the  momentum  space  density  matrix  elements,  pie(p,  p7;  t)  =  ple(p,  p7; 
by 


p'u(R,t)=p[ee-‘^R-^‘1  (35) 

=  J  J  <iWrie(p.  P';  p  -  p'  +  fikp).  (36) 

The  coherence,  p'le(p,  p7;  t)einpt,  has  been  written  in  the  Schrodinger  representation  and  is 
obtained  from  the  density  matrix  for  the  atom  plus  pump  fields  by  tracing  over  the  number 
of  photons  in  the  pump  fields, 


Pie(P’  p';  t)  =  e  l^pte  luJlet  ^2  Pie(p>  n2!  p7,  ni,  n2;  t).  (37) 

ni,ri2 

where  p[e( p,  ni,  n2\  p7,  t)  is  in  the  interaction  representation  with  respect  to  the  inter¬ 

nal  energy  levels  and  pump  fields.  One  cannot  derive  a  differential  equation  for  ple(p,  p7;  t ) 
starting  from  the  original  Hamiltonian  in  Eq.  (I3aj)  since  the  excited  state  has  been  adiabat- 
ically  eliminated  from  the  effective  Hamiltonian.  However,  by  reintroducing  an  interaction 
term 


Haf  =  h 


gi  \e)  (1|  e*kl'R  +  g2  \e)  (2|  e‘k2H  +  Xp  \e)  (1|  e 


ik2-R 


si(ki-R-f2pt) 


+  h.c. 
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into  the  Hamiltonian  and  writing 


Pie  (P?  p';  t)  =  (  e  *Qpfe  2a;iet  p[e(  p,  m,  ra2;  p',  m,  n2;  t)  j  .  (38) 

V  ni,n2  / 

one  finds  the  equation  of  motion  for  p[e( p,  ni,  n2;  ph  n( ,  rtf,)  to  be 

Pie(P)  wi,  n2;  p,  n),  n'2)  =  H<nPp/  -  (7  +  7e)/2]  pL(P,  w2;  p',  ni,  n'2) 

-ix*ie<Q^ei)t  [p[e(p  +  nkl5  m  -  1,  n2;  ph  ni,  n2) 

~Pn(p,  wi,  n2;  p'  -  hkun\  +  1,  n'2)\ 

-»Xpet(np-Wel)t  [p^e  ( P  +  ^kP>  ™i,  n2;  p,  n( ,  n2) 
-p{i(p,ni,n2;p/  -  Pkp,  n(,  n'2)] 

+*X2e*(°2+a;2e)Vi2 (P?  ni,  n2; p'  -  hk2,  ni,  n'2  +  1),  (39) 


where  the  t  argument  has  been  dropped.  By  carrying  out  the  trace  in  Eq.  (371),  one  obtains 
terms  such  as, 

Pn(p:P  -  frki)  =  Pii(P>ni»n2;p/  -  hki,m  +  1,  n2), 

ni,ri2 

so  that  the  equation  of  motion  for  p/le(p;p/)  has  a  form  which  is  identical  to  that  which 
would  have  been  obtained  using  classical  pump  fields, 


Ple(p;  p)  =  -  [*(fip  +  Wle  +  WPP')  +  (7  +  7e)/2]  pie(P;  P) 
-IXp  [Pee(  P  +  frkp!  p)  -  Pn(p;  P  -  fckp)] 

[p^e(p  +  hkp,  p')  -  pn(  p;  p  -  ftki)] 

+*X2e_*<5p2Vi2(P!  p'  -  ftk2). 


In  terms  of  the  perturbation  series  solution 


Pn(p;p') 

Pi  2(p;p') 


Pn  (p;  p)  +  Pn(p;  p)e^lf  +  Pn(P:  p')^1*; 

pS?(p;  p)  +  Pi2(p;  p'K5plt  +  pr2(p;  p')e_l5plt 


e 


i(5i2i . 
> 


(40) 


(41a) 

(41b) 
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where  p^(p;p')  are  independent  of  Xp  and  pfj>{ p;  p  )  are  linear  in  Xpi  the  steady  state 
solution  for  large  detuning  is 


p')  « }  (x>£(p; p'  -  nk2)  +  U/>n ip:  p'  -  fiki)  +  x’pPi. ?(p; p' 


(42) 


By  making  a  change  of  variables,  p  —  p;  — s ►  p  —  p'  —  Kkp  in  Eq. 


one  gets, 


Ple(R^)  = 


g  i(kp-R 

(2nh)3A 


d3pd3p'ei{p-p')n/h5( 


P-P 


x{x 2  (Pr2  (P-  p'  +  hkp2)  +  pf2{p-hkpl;  p'  +  hkia 


0)  + 


( in-  lV  4-  fik 


(43) 


One  must  now  obtain  equations  for  p(1(p;  p')  and  p{2(p;  pO  using  the  effective  Hamiltonian 
Htot=  H  +  V/,  solve  these  equations,  and  then  extract  p^p;  pO  and  Pn(p;  pO  from  these 
solutions  using  Eqs.  (I41a;-I41b{) . 

The  terms  appearing  in  Eq.  (143;)  may  be  expressed  in  terms  of  the  dressed  state  density 
matrix  elements.  If  one  is  interested  only  in  terms  linear  in  yp,  one  can  expand  the  dressed 
state  density  matrix  to  first  order  in  Xp  as 


Pir 


Pi  v  f  Pi i' 


(44a) 


Pjj' 


10420* . 


(44b) 


where  /,  J  =  {H0,H0},  I'  =  {Hi,  Hi},  and  J'  =  {H2,H2}.  The  pd.( p;  p')  needed  in  Eq.(l43|) 
can  be  expressed  in  terms  of  dressed  state  density  matrix  elements  using  Eqs.  (211,221,149!)  as 
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Pn(p;  P  +  ^kpi)  =  sin  00  (p)  sin0i(p)p^oAi 


sin  6*o (p)  cos  0i(p )Pa0b1 


(45a) 


-  cos  0O (p)  sin  0i(p)PsoAl  +  cos  0o(p)  cos  91  (p )PbqBi  5 
Pn(p-^kpi;  p)  =  sin 6*2 (p)  sin0o(p)PA2Ao  “  sin02(p)  cos  0o(p)pXbo 

—  cos  6*2 (p)  sin  9o(p) Pb2a0  +  cos  02(p)  cos  0o(p)Pb2.bo;  (45b) 

Pn( P;  P  +  ^2)  =  el(pd  (-  sin  0o(p)  cos  0i(p)/4oa1  -  sin0o(p)  sin  0i(p)/oJ,Bl 

+  cos0o(p)cos0i(p)p£oAi  +cos0o(p)sin0i(p)p^oBi)  ;  (45c) 

Pi2(p-^kpi;p  +  ^ki2)  =  (-sin02(p)cos0o(p)PA2Ao  -  sin02(p)  sin0o(p)PA2So 

+  cos  6*2 (p)  cos  9o(p) Pb2Aq  +  cos02(p)  sin0o(p )Pb2b0)  •  (45d) 

The  state  vector  in  the  Schrodinger  representation  may  be  expanded  in  terms  of  the 
dressed  states  for  the  0,  1,  and  2  manifolds, 


In  the  following,  the  momentum  labels  are  suppressed.  The  Schrodinger  equation  for  the 
dressed  state  amplitudes  is  then  given  by 


ih c  =  (H0  +  Vm)c 


(46) 


where 


(47) 


and  matrix  elements  of  V [rj  are  given  in  Eqs.  (33aJ-l33hj) . 


Using  Eq.  (46;)  along  with  the  ground  state  decay  rate  7  and  incoherent  pumping  to  the 
0  manifold,  one  finds  that  density  matrix  elements  for  the  six  dressed  states  in  the  three 
manifolds  (0),  (1),  (2),  evolve  as 


(48) 
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where 


Pd(p,p')  =  c(p)ct(p/) 


/ 


Pa,  a. 

Pa,b, 

Pb,a, 

Pb,b, 

Pa0a, 

Pa0b, 

Pb0a, 

Pb0b, 

Pa2a, 

Pa2b 1 

y  pb2a. 

Pb2b, 

Pa,a0  Pa-,  b0 
Pb,a0  Pb,b0 
PAqAq  Pa0bo 
Pb0a0  Pb0bo 
PA2Aq  PA2B0 
Pb2a0  Pb2Bo 


and  the  pumping  matrix,  A d,  has  the  block  diagonal  form 

Ad  =  diag(0,  Ad(p,  p;),  0), 


Pa,a2 

Pb,a2 

Pa0a2 

Pb0a2 

Pa2a2 

Pb2a2 


Pa,b2 

Pb,b2 

PAqB2 

Pb0b2 
Pa2b2 
pb2b2  y 


(49) 


(50) 


where  Ad(p)  has  the  basic  structure  given  in  Eq.  (l24j),  modified  to  allow  for  p,p7  coher¬ 
ence.  In  particular,  the  off  diagonal  elements  of  Ad(p,  p7)  that  give  rise  to  the  nonsecular 
contribution  to  the  line  shape  are  of  the  form 


Aab(p,  p'+hku)  =A2(p+hki2,  p'+hk12)  cos  [0o(p)]  sin  [6>o(p7)]-Ai(p,  p')  sin  [0o(p)]  cos  [0o(p')] . 


The  general  form  of  the  solution  is  linked  to  the  incoherent  pumping  of  levels  1  and 
2.  For  a  subrecoil  cooled  vapor,  the  pumping  rate  density  for  bare  state  density  matrix 
elements  p,j( p;  p7)  is  assumed  to  be 

Aij(p,  p')  =  AjV-1(27i7i)3<5(p)<5(p  -  p7)^-  (51) 

where  A*!/-1  has  the  dimensions  of  (volume x time)-1  and  can  be  interpreted  as  the  pumping 
rate  to  state  1  or  2  in  position  space.  With  this  form  of  pumping,  ph(p:  p7),  must  be 
proportional  to  <5(p  —  p7  —  Kk')  where  k7  is  some  algebraic  combination  of  the  pump  and 
probe  field  propagation  vectors.  To  obtain  p'le  from  Eq.  (j43|),  one  must  keep  only  those 
terms  in  the  integrand  of  Eq.  (|43|)  proportional  to  <5  (p  —  p7)  .  The  pg}(p;  p)  term  in  Eq. 
(•43;)  makes  no  contribution  to  the  absorption  since  it  is  real  and  will  be  ignored  from  this 
point  on. 
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The  incoherent  pumping  of  states  1  and  2  populate  two  different  manifolds.  The  pumping 
of  state  1  populates  the  (p  =  0,rii,ra2)  manifold  while  the  pumping  of  state  2  populates  the 
(p  =  — hki2,  n±,  n2)  manifold  since  this  manifold  involves  state  2  with  zero  momentum.  Thus 
in  viewing  absorption  or  emission,  two  distinct  initial  state  manifolds  must  be  included, 
leading  to  the  possibility  of  sixteen  rather  than  eight  components  of  the  spectrum.  Here  we 
set  A2  =  0. 

Substituting  Eqs.  (|45a|-f45d;)  into  Eq.  (143;)  and  using  the  steady  state  solutions  for  pfp 
and  pjj,  obtained  from  Eqs.  (!44a|-|44b{)  and  Eq.  (148;),  one  finds,  after  some  manipulation, 
the  final  expression  for  the  absorption  coefficient, 


a 

«o 


= 


a 

«0  ,  sec 
\G\Aab 


+ 


a 

a0y  ns 

sin  26*o 


-  (  i>Vs in  0i  -  -  cos  6*i  )  (TaqA^A')  +  rSoAl(A')) 


V«o  )ns  I A  |  7  2((4°i>)2  +  72) 

-  cos  6>1  +  i  sin  (r^s^A')  +  TSoSl(A')) 

+  p2  -  r2  -  2^  cot  2$o )  (sin2  82  (rJoj2(A')  +  rBo242(A'))  +  cos2  ft,  (rW  A')  +  rSos2(A'))) 


rWA')  = 
rBoJl(A')  = 
CWA')  = 

TsoSifA')  = 

rWA')  = 

rBoyi2(A')  = 

r^B2(A')  = 
rSos2(A')  = 


A' 


1 


LUhk. 


7^ AB  ~  U AB )  )  ^AB  +  f  LA0A1  (A'); 


,(1) 


-pi 


,(0) 


A' 


LUhk. 


.<°>  u-  o;(1) 
,(i) 


A'  -  UKkpl  +  -  (UaB  +  u AB )  )  ^AB  +  72 


-pi 


1  (, ,(°)  4- ,  ,w  \  \ ,  1  y 

2  y^AB  +  u AB  )  )  u AB  +  7 


1  A'  1  /  (0)  (1)  \  \  (Uj  .  2 

-  (  A  -  LUjik.pi  +  2  \UAB  ~  UAB>  ]  u AB  +  7 


-  (  A'  +  UJjikpi  +  ^  (^AB  _  w AB )  )  ^AB  +  72 


;(1) 

'AB 

,(2) 


,(°) 

;AB 

,(0) 


A'  +  LUh k  1  -  r(^AB  +  uab)  )  02 AB  +  7 


Kpi  2 


(2) 
AB 

(°)  ,  ,  ,(2) 

AB 

(0)  ,  (2) 


,(0) 


-  (  A'  +  tu^kpi  +  2  (UAB  +  UAb)  )  WAB  +  72 


A'  +  ^kpl^-Ki-^))^+7 


,(o) 


LBoAi(A'); 
LaqB^A'); 

LboBi(A'); 

La0a2{ A  ); 
Lb0a2(A/); 

La0b2(A'); 
Lb0b2(A/); 


where  Aab  =  —  |Ai  sin  (20o)  and  (  —  )  is  given  by  Eq.  (|2?aj).  The  reason  the  absorption 

V  a°  /  sec 

coefficient  is  expressible  as  a  sum  of  the  secular  term  plus  a  non-secular  term  is  linked 
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to  the  fact  that  the  secular  approximation  consists  solely  of  neglecting  the  off-diagonal 
components  of  A 4.  Since  the  first  order  solutions,  p\v  and  Pjj/,  are  linear  in  the  pumping 
terms,  (  —  )  contain  terms  proportional  to  and  A b  while  (  —  )  is  proportional  to 


A ab-  This  simplification  would  not  occur  for  a  more  complex  decay  scheme  for  states  |1  > 
and  1 2  >  since  the  decay  would  couple  density  matrix  elements  in  a  field  dependent  manner 
(see  [g3]). 

The  nonsecular  term,  (  )  ,  consist  of  dispersion-like  structures  centered  at  the  same 

V  a°  /  ns 


locations  as  the  resonances  in  (  — 

<*o 


.  In  the  secular  limit, 


QO 


•C 


OLQ 


and  ( -a 

a0 


can  usually  be  ignored.  Notice  that  if  one  chooses  ij>  and  rj  such  that  a  pair  of  absorption 


lines  in 


ao 


vanish,  then  the  corresponding  terms  in 


a0 


also  vanish  so  that 


is 


identically  zero.  However,  this  will  not  be  true  for  the  gain  terms  in  (  ^ 


a°^ 

since  the  cor- 


sec 

responding  terms  in  (  —  )  have  a  different  interference  coefficient,  {rj2  —  r]~2  —  2^  cot  2$o)- 


\  /ns 

Figure  6  shows  a  plot  of  the  non-secular  absorption  coefficient  for  the  same  parameters  as 
Fig.  3.  In  this  plot,  the  non-secular  terms  are  ~  1000  times  smaller  than  the  secular  terms. 
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Figure  1.  Schematic  diagram  of  atom-field  system. 

Figure  2.  (a)  Transitions  between  the  0  manifold  and  the  1  and  2  manifolds  leading  to 
probe  gain  or  absorption  in  the  bare  state  picture.  The  states  in  the  1  and  2  manifolds 
are  displaced  from  the  states  in  the  0  manifold  by  an  amount  ~  MR.  The  Rabi  frequencies 
shown  are  those  which  couple  the  states  from  the  0  manifold  to  the  1  and  2  manifolds,  (b) 
Energy  levels  in  the  dressed  state  basis.  The  center  of  the  1  manifold  has  an  energy  tko\o 


above  the  center  of  the  0  manifold  and  similarly,  the  center  of  the  2  manifold  is  tk02o  below 
the  0  manifold,  (c)  Illustration  of  the  coupling  of  the  pump  and  probe  fields  to  the  dressed 
states  for  the  |Ao  >— >  \A\  >  transition.  The  Rabi  frequencies  shown  are  those  that  couple 
the  bare  state  |1  >  component  of  |Ao  >  to  the  bare  state  components  of  \A\  >. 

Figure  3.  Plot  of  (  —  )  for  A2  =  0,  A1/7  =  1,  ip  =  —1  and  77  =  2  showing  all  eight 
absorption  and  emission  lines.  The  detuning  is  R2/7  =  300  and  the  two-photon  pump 
Rabi  frequency  is  \G\/ry  =  250.  The  recoil  energies  are  ^nkl2 !l  —  40,  07^/7  =  60,  and 
frkpi  ■  ki2/M7  =  80.  Note  that  7  -C  ye  where  ye  is  the  excited  state  decay  rate. 


Figure  4.  Plot  of 


showing  destructive  interference.  Parameters  are  the  same  as 


Fig.  4  except  for  ip  and  77.  The  solid  line  corresponds  to  ip  —  —1  and  77  =  \J tan  6\  =  0.8383 
while  the  dotted  line  corresponds  to  ip  —  +1  and  77  =  \J cot  6\  =  1.1928. 

Figure  5.  Plot  showing  sixteen  lines  in  the  (  —  )  probe  spectrum  for  A2/7  =  A1/7  =  1, 
ip  =  + 1  and  77  =  0.1.  The  detuning  is  $12/7  =  500  and  the  two-photon  pump  Rabi  frequency 
is  \G\/ry  =  750.  The  recoil  energies  are  =  50,  07^1/7  =  75,  and  hkpi-k^/My  =  —50. 


Figure  6.  Plot  of  for  the  same  parameters  as  Fig.  (3).  The  non-secular  absorption 

coefficient  has  dispersionlike  structures  at  the  same  location  is  the  line  centers  of 
The  amplitudes  of  these  non-secular  terms  is  typically  ~  1000  times  smaller  than  the  secular 
line  strengths,  consistent  with  7 /^ab  =  0.00177. 
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